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involved. For Alfven waves to propagate a^ong cached field lines, 
the Binimum rest nass density becomes p > c ) = 40 g cm , a 

value that violates the assumption of a vacuum magnetic field. 
Similarly if one postulates a high temperature region free of mag- 
netic field as would be required by Fenimore et al 7 to Compton heat 
a soft photon flux, one has the problem of how to keep these 
electrons hot using a source of energy that does cot introduce more 
electrons than is allowed by the "thin" approximation. 

For these reasons we have spent some time attempting to find an 
emission mechanism for gamma bursts that converts a reasonable 
source of energy into the hard photon spectrum observed. We believe 
we have found such a mechanism as the normal consequence of moder- 
ately thin (l = 1 to 4) accretion onto a field free neutron star. 
The model is that photons of nearly constant number are sequentially 
compressed, heated, diffused, further compressed, heated, diffused 
and ultimately escape in a hard spectrum at an average flux near the 
Eddington limit, but with wide fluctuations. 

The physics of such accretion is complex. To begin with: 

1. The mass flux corresponding to steady state Eddington limit 
accretion corresponds to a thickness of 1=2. The accreting 
matter cannot be much thicker, else the energy flux of radia- 
tion will be less than the kinetic energy flux of infalling 
matter . 

2. In a transient state when the initial instantaneous photon flux 
corresponds to a value much less than the Eddington limit flux, 
<|>Edd* the free fall time is very much less than the time neces- 
sary to multiply the photon number by bremsstrahlung. There- 
fore, accretion with I = 2 to 4 occurs without increase in the 
photon number. 


CALCULATIONS 

We have modeled this phenomenon with a one-dimensional (radial) 
multi-energy group numerical calculation with mass layers and photon 
diffusion in both energy ai_d radius. Both radiation stress and 
gravity accelerate the accreting matter. The neutron star surface 
is assumed to have an albedo of unity, mass layers are merged with 
the surface when they get within AR of it and each layer is treated 
as a semi-transparent diffuse scattering layer with zero heat 
capacity and therefore no defined temperature. A photon number 
density corresponding to a small temper ire initially seeds the 
problem. A typical resulting spectrum is shown in the figure where 
the initial seed temperature was 0.355 keV, ^ = 2.2 (integrated 
to infinity) and the initial velocity distrit?uVion corresponded to 
free fall. The curve drawn through the points is that corresponding 
to a thin bremsstrahlung distribution with kT g = 139 keV. 


involved. For Alfven waves to propagate a^.ong cached field lines, 
the ninimum rest mass density becomes p > c ) = 40 g cm , a 

value that violates the assumption of a vacuum magnetic field. 
Similarly if one postulates a high temperature region free of mag- 
netic field as would be required by Fenimore et al 7 to Compton heat 
a soft photon flux, one has the problem of how to keep these 
electrons hot using a source of energy that does not introduce more 
electrons than is allowed by the "thin" approximation. 

For these reasons we have spent some time attempting to find an 
emission mechanism for gamma bursts that converts a reasonable 
source of energy into the hard photon spectrum observed. We believe 
we have found such a mechanism as the normal consequence of moder- 
ately thin (t = 1 to 4) accretion onto a field free neutron star. 
The model is that photons of nearly constant number are sequentially 
compressed, heated, diffused, further compressed, heated, diffused 
and ultimately escape in a hard spectrum at an average flux near the 
Eddington limit, but with wide fluctuations. 

The physics of such accretion is complex. To begin with: 

1. The mass flux corresponding to steady state Eddington limit 
accretion corresponds to a thickness of 1 = 2. The accreting 
matter cannot be much thicker, else the energy flux of radia- 
tion will be less than the kinetic energy flux of infalling 
matter . 

2. In a transient state when the initial instantaneous photon flux 
corresponds to a value much less than the Eddington limit flux, 
4>Edd» t ^ e free fall time is very much less than the time neces- 
sary to multiply the photon number by bremsstrahlung. There- 
fore, accretion with I = 2 to 4 occurs without increase in the 
photon number. 


CALCULATIONS 

We have modeled this phenomenon with a one-dimensional (radial) 
multi-energy group numerical calculation with mass layers and photon 
diffusion in both energy ai_d radius. Both radiation stress and 
gravity accelerate the accreting matter. The neutron star surface 
is assumed to have an albedo of unity, mass layers are merged with 
the surface when they get within AR of it and each layer is treated 
as a semi-transparent diffuse scattering layer with zero heat 
capacity and therefore no defined temperature. A photon number 
density corresponding to a small temper ire initially seeds the 
problem. A typical resulting spectrum is shown in the figure where 
the initial seed temperature was 0.355 keV, ^ = 2.2 (integrated 
to infinity) and the initial velocity distribution corresponded to 
free fall. The curve drawn through the points is that corresponding 
to a thin bremsstrahlung distribution with kT g = 139 keV. 



log energy (key) 


Fig. A Y“burst spectrum from field free matter accreting 

onto a neutron star. The curve is a bremsstrahlung fit 

with kT = 139 keV. 
e 

The justification for the set of assumptions leading to the 
calculations are: 

1. In steady state Eddington flux limited accretion the gravita- 

tional field of the neutron star acting on the ions drags the 
electrons by a ^ charge separation electric field E < g(Jm /e 
where g = (J is the molecular weight per electron, ind 

m the mass of the proton. This electric field, E = 100 v/cm 
ii small enougji^ that charge neutrality almost prevails, 

An /n = 3 x 10 . This ensures that the electrons are 

0 0 

dragged by the ions against the photon stress even for stopping 
distances that are infinitesimal compared to the neutron star 
radius . 

2. The seed photon number density, generally characterized by a 

black body temperature kT. . < 1 keV is small compared to the 

Eddington limit value of = 2 keV for a neutron star of mass 
1 M_. This ensures initial free fall of the plasma. 
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The photon number density always greatl^exceeds the e^ctron 
number density since n is of„j>rdec 1(T i/R^ or 10* cm 
while at 1 keV, n^ y = f.5 x 10 cm Therefore the electron 
specific heat is negligible compared to that of the photons. 
Provided the electrons do not run away in the charge separation 
electric field, the electron temperature will be an average of 
the local photon energy. Hence electron recoil and comptoniza- 
tion where the thermal energy content of the hot or cold elec- 
trons heats or cools the photons can be neglected In fact, 
runaway will not occur because an electron undergoes many 
photon collisions in a distance equal to its energy multiplica- 
tion distance \ - kT /Ee, Ee < n. < hv> o /3t, t = 

2A/(2 kT /m ) , and hence tne number of collisions in C \ becomes 
n hv °c C ^. ~ e (c/v e ) sufficient to ensure against runaway. This 
assumes that kT = <h\»/3 consistent with the large number of 
collisions . 

Since for matter freely falling at a constant rate, /T pdR « 
R , the thickness pR ^ R ^ and so photons will scatter and 
be heated within several radii of the neutron star. From this 
distance the free fall time is roughly equal to the diffusion 
time. This ensures that some photons are trapped and raised in 
energy by compression, yet also ensures that leakage allows the 
large photon heating without exceeding the ^ddington limit. 
Since the gravitational potential is roughly c /5, the limiting 
free fall velocity is = c/2. The diffusion "velocity" is 
roughly (c/3l) so that t = 1 to 2 meets the criterion. 

The neutron star surface has an albedo of the order of unity. 
Since bremsstrahlung is small compared to Compton scattering, 
the "surface" is where there first exists an additional scat- 


tering mean free path. Therefore a photon reaching the 
"surface" from larger radius has a high probability of being 
scattered back out, rather than being absorbed and reemitted as 
is necessary to establish a black body spectrum. 

The photons added to the accreting matter during free fall due 
to the neutron star surface emission are not large compared to 
the assumed in situ number density. Since the free fall time 
is comparable to the escape time, the net addition of photons 
after the start of free fall is net greater than several fold. 
Photons that remain local, i.e., trapped, are compressed by a 
factor no greater than the local matter compression ratio, r|, 
so that photons that remain trapped within a -Aiyen layer of 
matter are increased in energy by the factor r| 7 . This is a 
constraint in the numerical calculation. If the problem con- 
sists of n layers and a minimum R/AR = rj » the maximum 
increase, in energy of a given photon becomes of the order 
(R/AR) n . 

The Eddington limit flux just balances radiatiap pressure 
against the gravitational force, which is GM^«,/R per unit 
mass. Tc decelerate freely falling matter in a distance AR 


requires an impulse per unit mass of the momentum v be 
supplied in a time 2AR/v and therefore requires a force v /2AR. 
The ratio of the flux required to decelerate the matter to the 


2 2 

Eddington flux is therefore v R /2GMj <s AR which is R/AR taking 
account of v 2 = 2GM^/R. The bremsstrahlung rate is propor- 
tional to n 2 or AR and the time of deceleration is propor- 
tional to AR so that the total number of photons produced by 


bremsstrahlung is proportional to 1/AR the same as the decel- 
eration flux. Therefore, for all AR and the previous con- 
ditions, the photon heating rate by the free fal] matter will 
exceed the photon multiplication rate by bremsstrahlung. 

10. Double corapton scattering will increase the photon number at a 
rate compared to Compton scattering: 
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The highest energy photons have on the average undergone the 
most scatterings and so are the most likely to have made a 
second photon along the way. The number of scatterings per 
doubling is- roughly 3, so the number of scatterings going from 
(3 keV -*■ me 2 ) = 3 log- (160) = 22. Then 
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which is < 1 so that double Compton does not limit photon 
energy significantly until h\> > me . 

11. The free falling matter is decelerated only when the photon 
energy exceeds the matter energy. Neither bremsstrahlung nor 
double Compton scattering affect the total photon energy, the 
first because of limited time, and the second because it merely 
redistributes the photon energy, i.e., mal.es a softer spectrum. 
We have chosen our thickness such that compression heating 
proceeds at just about the rate that photon energy is lost by 
diffusion. To illustrate this point we take the diffusion 
photon density distribution between at a reflector at x = 0 and 
an outer boundary at X: 

p = cos (tix/2X) 

The collision rate is 

N coll = Jo COp dx 


and the escape flux is 
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The escape per collision is n /(12i ). Hence if it requires 3 
collisions on the average to double the photon energy, then the 
total retained photon energy will remaia. constant if the prob- 
ability of escape per collision is 2 , corresponding to 
1=1. Hence a modest increase in I, will effectively trap the 
photons and their .Jtotal energy can increase without limit 
proportional to q . Hence we feel a natural limit to com- 
pression will be reached when collisions lead to pairs that 
then substantially increase I. This offers a natural photon 
energy limit to the spectrum that is flux sensitive. 

It seems reasonable to assume that the photon diffusion 
terminates when pair production doubles the lepton density. 
The calculations do not include pair production or this limit. 

CALCULATION PARAMETERS 

The calculation of the figure is relatively _insensitive (= 20%) 
to changes in energy group size (W /W.) < Jl, number of zones 
n > 10, and AR/R < 0.005. If the initial 1 seed temperature of .35 
keV is reduced, then the energy group size, and AR must be smaller 
and n greater to obtain a result not grossly sensitive to the dis- 
creteness of variables. 


ASTROPHYSICS 

The maximum magnetic field that can be pe^.-itted can be esti- 
mated from Liang 5 who estimates the electron density for equal 
emission from bremsstrahlung and synchrotron radiation to be 

n g = 2.3 x 10 31 T 3/2 F(T) Bj 2 Z _1 cm" 3 , 

2 

where T is kT /me and F(T) is of order unity for T = 1. Hence to 
keep synchrotron, radiation below bremsstrahlung in the low density 
plasma (n p = 10 cm 3 ) requires <10 , or B < 3 x 10 gauss. 
Since bremsstrahlung is rather a weak limit on the process and the 
high energy photons come from compressed layers, a field of several 
10 gauss would be acceptable. However, the limiting ^.eld is far 
less than what would control or funnel accretion, 3 x 10 gauss. In 
any event magnetic field precludes the usual interpretation of the 
spectral structure as due to a cyclotron line. Some of our calcu- 
lations show features similar to the observed dips at a few tens of 
keV. Another possibility is that a parallel softer photon source 
gives rise to these features. 11 

We feel that one possible astrophysical circumstance in which 
the field-free accretion we describe could occur is the instability 
of a Keplerian disc when it intersects, i.e. touches, the neutron 
star surface. Friction and turbulence will then cause some low 


density nass to accrete. Other possible candidates are super- 
Eddington thermal conduction to the surface from the interior of the 
neutron star or shocks that sequentially blow off surface layers 
that subsequently fall back. In both these possible examples one 
difficulty is the high atomic number of the neutron star matter 
results in high opacity and high photon emission. Comet accretion 
into an orbit that subsequently decays by gravitational wave 
emission is another possibility. 
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